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ABSTRACT 


The  bulk  second  order  closure  method  of  Garwood 
(1977)  is  used  to  predict  the  well-mixed  layer  depth  and 
temperature  evolution  at  Ocean  Station  Papa  for  the  years 
1953-1969.  The  time  dependent  boundary  conditions  are  the 
specification  of  the  surface  buoyancy  flux  and  wind  stress 
derived  from  the  three-hourly  atmospheric  observations.  Much 
cf  the  variance  of  the  well-mixed  layer  depth  and 
temperature  evolution  is  obviously  associated  with  the 
annual  cycle.  However,  a  closer  inspection  of  the  results 
reveals  a  variability  related  to  the  synoptic  response  of 
sea-surface  temperature  and  mixed  layer  depth  to  both  strong 
and  liqht  wind  events.  The  simulations  are  not  intended  to 
be  a  best  fit,  but  rather  a  demonstration  of  the  capability 
of  a  one-dimensional  model  to  simulate  the  interannual 
variability  observed  at  OHS  papa. 


2 


1.  INTRODUCTION 


Ocean  Station  Papa  (50N,14  5W)  in  the  eastern  North 
Pacific  Ocean  has  been  the  site  of  a  number  of  one- 
dimensional  mcdel  simulations  of  the  oceanic  planetary 
boundary  layer  or  mixed  layer.  Using  the  bulk  model  of  Kraus 
and  Turner  (1967)  at  this  site,  Denman  and  Miyake  (1973) 
first  simulated  the  mixed  layer  response  to  local 
atmospheric  forcing  during  June  1970  for  a  twelve-day 
period.  Mellcr  and  Durbin  (1975)  applied  a  profile  (rather 
than  bulk)  turbulence  closure  model  to  simulate  a  few  weeks 
of  upper  ocean  thermal  structure  evolution  for  the  same  time 
of  the  year  (late  spring  but  after  the  mixed  layer  had 
shallowed).  Camp  and  Elsberry  (1978)  used  three  versions  of 
the  bulk  mixed  layer  model  (Elsberry  et  al ,  1976;  Kim,  1976; 
and  Kraus  and  Turner,  1967)  to  study  periods  of  fall  and 
early  winter  deepening  in  response  to  strong  atmospheric 
forcing  at  a  number  of  locations,  including  Ocean  Station 
Papa.  The  lengths  of  these  model  integrations  were  up  to  a 
month,  and  a  number  of  different  years  were  examined.  The 
Mixed  Layer  Experiment  (MILE)  of  August  18  to  September  5, 
1977  was  also  situated  at  Ocean  Staion  Papa.  At  least  two 
different  model  simulations  have  bee;  conducted  of  the  MILE 
period;  Garwood  (1978)  and  Davis  et  al  (1982).  All  of  the 
above  model  simulations  at  Ocean  Station  Papa  were  of 
relatively  short  duration,  and  all  were  during  the  same  half 
of  the  year,  namely  between  the  months  of  May  and  December 
after  the  seasonal  thermocline  had  been  set  up. 


For  abcut  30  years,  from  the  early  1950's  until 
1981,  the  Canadian  Goverment  sponsored  a  nearly  continuous 
collection  of  meteorological  and  oceanographic  data  at  Ocean 
Station  Papa.  The  Canadian  Goverment  has  recently 

terminated  ship  activities  at  Papa,  thus  interupting  one  of 
the  few  data  sets  available  to  oceanographers  who  wish  to 
study  lonq  term  variability  in  the  upper  ocean.  Such 
observations  have  included  BT  drops,  air  temperature , 
humidity,  wind  measurements  and  shy  cover. 

The  purpose  of  this  study  is  to  report  the  multi¬ 
year  evolution  of  the  ocean  thermal  structure  predicted  by  a 
boundary  layer  entrainment  model  in  response  to  forcing 
parameters  calculated  from  observations  taken  at  OIS  Papa. 
In  spite  of  the  existence  of  this  long  time  series  of 
observations,  until  now  no  model  simulations  have  been 
reported  which  are  for  more  than  a  few  weeks  in  duration, 
and  no  study  has  previously  been  conducted  of  the 
interannual  variability  in  the  ocean  mixed  layer  response  to 
the  atmosphere.  Also,  no  simulations  have  previously  beer, 
conducted  for  the  period  of  spring  transition  (usually 
between  Harch  and  Hay  at  Ocean  Station  Papa)  when  the 
seasonal  thermccline  is  set  up.  The  purpose  of  this  paper  is 
to  document  the  results  of  a  series  of  seventeen  consecutive 
ncdel  integrations  (each  of  one  year  in  duration)  that  for 
the  first  tine  take  full  advantage  of  the  continuity  of  the 
data  set. 
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2.  METHOD 
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The  bulk  second  order  closure  method  of  Garwood 
(1977)  is  employed  as  the  OPEL  model.  The  required  time 
dependent  surface  boundary  conditions  is  the  spec! f icat f on 


of  the  surface 

buoyancy  and 

momentum  flux 

which 

are 

calculated  using 

bulk  aerodynamic  formulae 

(Camp 

and 

Elsberry, 1978) . 

The  forcing 

is  calculated 

from 

the 

observations  taken  every  three 

hours.  Except 

for 

the 

parameterization  of  the  shear  production  of  turbulence, 
there  is  r.o  calculation  of  the  momentum  budget  in  this 
study.  For  each  c*  the  seventeen  years,  the  model  is 
initialized  with  a  mixed  layer  depth,  mixed  layer 
temperature,  a  temperature  'jump  at  the  base  of  the  mixed 
layer  and  a  bottom  temperature  (in  this  case  200  m)  tc  which 
the  temperature  decreases  linearly  from  that  temperature 
lust  below  the  initial  mixed  layer  depth  (see  Table  I  for 
these  values).  Each  cf  the  model  runs  is  also  initialized 
with  a  salinity  profile  that  is  a  constant  34  ppt  down  to 
the  mixed  layer  depth,  decreases  linearly  to  35  ppt  down  to 
a  depth  cf  150  m  and  remains  a  constnt  35  ppt  below  150  m. 
This  salinity  profile  contributes  to  the  late  winter  mixed 
layer,  but  has  no  influence  upon  the  evolution  of  the 
seasonal  thermccline  because  the  surface  salinity  flux  is 
zero.  This  is  equivalent  to  precipitation  balancing 
evaporation. 

Within  the  model  there  are  essentially  two  degrees 
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of  freedom  which  can  be  thought  of  as  “tuning  parameters" . 
For  all  seventeen  years  of  integration,  only  one  set  of 
tuninq  parameters  is  used.  These  tuning  parameters  were 
calibrated  in  such  a  way  so  that  the  maximum  summer 
temperature  predicted  by  the  model  for  1959  is  in  agreement 
with  the  maximum  summer  temperature  observed  at  OHS  Papa 
durinq  that  year.  The  year  1959  represents  no  special  cast 
tut  rather  a  random  choice.  This  was  not  intended  to  be  a 
best  fit  or  optimal  tuninq  in  this  study,  but  a 
demonstration  cf  the  relative  interannuai  variability  that 
can  be  simulated  by  a  cne  dimensional  model. 

3.  RESULTS 

Each  year  at  OHS  Papa  is  represented  by  a  set  of 
seven  fiqures  labeled  e-g.  Those  figures  marked  'a*  are  *he 
observed  mixed  layer  depths  from  BT  drops.  The  mixed  layer 
depth  was  defined  to  be  that  depth  where  the  water 
temperature  first  deviated  by  . 2  C  from  the  surface  value. 
Figures  marked  'b*  are  the  model  predicted  time  evolution  of 
water  temperature.  The  mixed  layer  depth  can  be  inferred 
from  the  tight  packing  of  isotherms.  Figures  marked  c-g  are 
a  four  day  running  mean  of  the  time  series  of  the  observed 
sea-surface  temperature,  air  temperature,  dew  point 
temperature,  wind  speed  and  percentage  of  cloudiness 
respectively. 

For  the  most  part,  the  seventeen  figures  are  self 
evident.  It  is  recccmmended  that  the  reader  peruse  them 


consecutively,  looking  for  interannual  differences  and 
similarities.  This  should  be  followed  by  a  much  closer 
inspection  of  each  year,  with  particular  attention  paid  to 
the  synoptic  response  of  sea-surface  temperature  and  mixed 
layer  depth  to  both  strong  and  ligh+  wind  events. 

Most  of  the  variances  in  the  dew  point,  sea  and  air 
temperatures  are  obviously  associated  with  the  annual  cycle. 
The  air  (dry  bulb)  temperature  in  all  of  these  cases  range 
from  winter  lows  of  about  4  C  or  5  C  to  summer  highs  ranging 
near  the  14  C  cr  15  C  mark.  On  the  average,  the  wind  speeds 
are  hiqher  in  the  winter  than  they  are  in  rhe  summer, 
although  the  highest  values  are  attained  in  the  late  fall. 
Oddly  enouqh,  the  percentage  of  cloudiness  is  at  a  minimum 
when  storm  activity  is  at  a  peak,  and  maximum  cloudiness 
occurs  in  the  summer. 

Since  the  model  predicted  values  were  only  saved 
cnce  a  day,  the  diurnal  signal  is  not  evident  here,  and  the 
annual  cycle  is  the  dominant  signal.  Initially,  the  mixed 
layers  are  deep  and  the  sea-surface  temperatures  are  low. 
Then,  sometime  in  sprinq,  the  upper-ocean  winter  regime 
qives  way  to  a  summer  regime  during  which  the  mixed  layer  is 
much  shallower  and  the  sea-surface  temperatures  are  much 
hiqher.  As  fall  approaches  and  ■•■he  frequency  of  storms 
increase,  the  mixed  layer  becomes  progressively  deeper  un  +  il 
once  aqain  the  winter  regime  is  restored.  There  is  a  limit 
to  mixed  layer  deepening  since  the  disspipation  length  scale 


prevents  the  mixed  layer  from  becoming  infinitely  deep.  The 
balocline  prescribed  as  an  initial  condition,  also  may  limit 
deep  mixinq  for  some  years,  as  revealed  by  model 
integrations  having  no  halocline  (not  shown  here). 

perhaps  the  most  striking  feature  in  the  year-to- 
year  variability  of  the  model  hindcasts  is  the  way  the 
winter  regime  is  transformed  into  a  summer  regime.  This 
phenomenon  is  sometimes  called  'spring  transition'.  During 
the  spring,  the  winds  are  lighter  and  there  is  an  increase 
in  the  ret  downward  heat  flux  into  the  mixed  layer.  The 
combination  of  the  lack  of  turbulent  kinetic  energy  to 
maintain  a  deep  mixed  layer  and  the  formation  of  a  layer  of 
buoyant  water  due  to  heating,  allows  a  new  mixed  layer  to 
reform  at  a  shallower  depth.  In  many  cases,  the  transition 
has  an  almost  sudden  occurrence  as  is  the  case  for  1959 
(Fig.  7a),  while  with  ether  cases,  the  change  from  a  winter 
to  a  summer  regime  is  more  gradual  and  occurs  in  two  or  more 
steps  as  in  1957  (Fig.  5a).  As  suggested  by  Elsberry  and 
Garwood  (1978),  the  date  of  spring  transition  strongly 
influences  summertime  sea-surface  temperatures.  An  early 
transition  means  that  heat  accumulates  in  a  shallower 
boundary  layer  for  a  long  period  of  time  and  hence  increases 
summer  sea-surface  temperatures  more  than  usual.  Conversely, 


a  late  transition  date  will  lead  to  relatively  cooler  summer 
temperatures.  At  OWS  Papa,  spring  transticn  normally  occurs 
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between  Julian  days  100  and  120. 

In  conclusion,  the  bulk  method  used  here  seems  to  be 
capable  of  simulating  the  seasoned  and  synoptic-scale 
variability  observed  at  OWS  Papa.  Seventeen  years  cf  air- 
sea  observations  and  model  integrations  provide  an  extensive 
test.  Certainly  other  effects  such  as  advection  would  be 
needed  to  provide  a  complete  explanation  of  ocean 
variability  at  OHS  Papa,  but  it  is  nevertheless  encouraging 
that  one  dimensional  mixing  can  account  for  a  large  part  of 
the  variance  on  time  scales  from  the  synoptic  to  more  than  a 
year. 
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TABLE  AND  FIGURE  CAPTIONS 


Table  I  Listing  of  initial  sea-surface  temperature,  mixed 
layer  depth,  temperature  jump  at  the  base  of  the  mixed  layer 
and  temperature  at  200  meters  for  the  years  1953-1969. 

Fig.  1.  (a)  Mixed  Layer  depths  from  BT  drops,  (b)  depth 
vs.  time  contours  of  temperature,  (c)  sea-surface 
teaperature,  (d)  air  temperature  (e)  dew  point  teapera+ure, 
ff)  wind  speed  and  (g)  cloud  cover  for  year  1953. 
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Table  I 


Depth 

Temp. 

T  jmp 

T  bot 

1953 

100.0 

6.6 

0.9 

WEM 

1954 

120.0 

5.4 

0.7 

3.4 

1955 

116.0 

6.9 

1.2 

3.4 

1956 

120.0 

5.5 

1.4 

3.0 

1957 

100.0 

6.7 

0.9 

mm 

1958 

120.0 

7.0 

0.9 

3.2 

1959 

106.0 

6.0 

0.6 
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1960 

110.0 

6.5 
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1.3 
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1965 
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5.0 

0.6 
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1966 
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1.0 

3.6 

1967 
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6.6 

1.4 

' 

3.2 

1968 

130.0 

5.3 

0.6 

3.2 

1969 
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5.6 

0.8 
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Figure  8.  (a)  Mixed  layer  depths  from  BT  drops, 

(b)  depth  vs.  time  contours  of  temperature, 

(c)  sea-surface  temperature,  (d)  air  temperature, 
(e)  dew  point  temperature,  (f)  wind  speed  and 
(g)  cloud  cover  for  year  1960. 


Figure  9.  (a)  Mixed  layer  depths  from  8T  drops, 

(b)  depth  vs.  time  contours  of  temperature, 

(c)  sea-surface  temperature,  (d)  air  temperature, 
(e)  dew  point  temperature,  (f)  wind  speed  and 
(g)  cloud  cover  for  year  1961 
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Rgure  10.  (a)  Mixed  layer  depths  from  BT  drops, 

(b)  depth  vs.  time  contours  of  temperature, 

(c)  sea-surface  temperature,  (d)  air  temperature, 
(e)  dew  point  temperature,  (f)  wind  speed  and 
(g)  cloud  cover  for  year  1962. 


Figure  It  (a)  Mixed  layer  depths  from  BT  drops, 

(b)  depth  vs.  time  contours  of  temperature, 

(c)  sea-surface  temperature,  (d)  air  temperature, 
(e)  dew  point  temperature,  (f)  wind  speed  and 
(g)  cloud  cover  for  year  1963. 
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Figure  12.  (a)  Mixed  layer  depths  from  BT  drops 

(b)  depth  vs.  time  contours  of  temperature, 

(c)  sea-surface  temperature,  (d)  air  temperature, 
(e)  dew  point  temperature,  (f)  wind  speed  and 
(g)  cloud  cover  for  year  1964. 
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Rgure  14.  (a)  Mixed  layer  depths  from  BT  drops, 

(b)  depth  vs.  time  contours  of  temperature, 

(c)  sea-surface  temperature,  (d)  air  temperature, 
(e)  dew  point  temperature,  (f)  wind  speed  and 
(g)  cloud  cover  for  year  1966. 
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Figure  17.  (a)  Mixed  layer  depths  from  BT  drops, 

(b)  depth  vs.  time  contours  of  temperature, 

(c)  sea-surface  temperature,  (d)  air  temperature, 
(e)  dew  point  temperature,  (f)  wind  speed  and 
(g)  cloud  cover  for  year  1969. 
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